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Abstract 

I review a variety of model systems and their quantum critical points, motivated by recent 
experimental and theoretical developments. These are used to present a general discussion of 
the non-zero temperature crossovers in the vicinity of a quantum critical point. Insights are 
drawn from the exact solutions of quantum critical transport obtained from the AdS/CFT 
correspondence. I conclude with a discussion of the role of quantum criticality in the phase 
diagram of the cuprate superconductors. 



Contributed chapter to the book 'Developments in Quantum Phase Transitions', 
edited by Lincoln Carr. 
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The author's book 1 on quantum phase transitions has an extensive discussion on the 
dynamic and transport properties of a variety of systems at non-zero temperatures above 
a zero temperature quantum critical point. The purpose of the present article is to briefly 
review some basic material, and to then update the earlier discussion with a focus on ex- 
perimental and theoretical developments in the decade since the book was written. We note 
other recent reviews^ from which portions of this article have been adapted. 

We will begin in Section UJ by introducing a variety of model systems and their quantum 
critical points; these are motivated by recent experimental and theoretical developments. 
We will use these systems to introduce basic ideas on the finite temperature crossovers near 
quantum critical points in Section [Til In Section [TTT| we will focus on the important quantum 
critical region and present a general discussion of its transport properties. An important 
recent development has been the complete exact solution of the dynamic and transport 
properties in the quantum critical region of a variety of (supersymmetric) model systems in 
two and higher dimensions: this will be described in Section IIVI The exact solutions are 
found to agree with the earlier general ideas discussed here in Section IIHI Quite remarkably, 
the exact solution proceeds via a mapping to the theory of black holes in one higher spatial 
dimension: we will only briefly mention this mapping here, and refer the reader to the article 
by Hartnoll in this book for more information. As has often been the case in the history 
of physics, the existence of a new class of solvable models leads to new and general insights 
which apply to a much wider class of systems, almost all of which are not exactly solvable. 
This has also been the case here, as we will review in Section |Vj a hydrodynamic theory 
of the low frequency transport properties has been developed, and has led to new relations 
between a variety of thermo-electric transport co-efficients. Finally, in Section I VI I we will 
turn to the cuprate high temperature superconductors, and present recent proposals on how 
ideas from the theory of quantum phase transition may help unravel the very complex phase 
diagram of these important materials. 

I. MODEL SYSTEMS AND THEIR CRITICAL THEORIES 
A. Coupled dimer antiferromagnets 

Some of the best studied examples of quantum phase transitions arise in insulators with 
unpaired S = 1/2 electronic spins residing on the sites, i, of a regular lattice. Using Sf 
(a = x, y, z) to represent the spin S = 1/2 operator on site i, the low energy spin excitations 
are described by the Heisenberg exchange Hamiltonian 
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FIG. 1: The coupled dimer antiferromagnet. The full red lines represent an exchange interaction 
J, while the dashed green lines have exchange J/X. The ellispes represent a singlet valence bond 
of spins (| U) - | lt»/V2. 

where > is the antiferromagnetic exchange interaction. We will begin with a simple 
realization of this model is illustrated in Fig. [TJ The 5 = 1/2 spins reside on the sites of a 
square lattice, and have nearest neighbor exchange equal to either J or J / A. Here A > 1 is 
a tuning parameter which induces a quantum phase transition in the ground state of this 
model. 

At A = 1, the model has full square lattice symmetry, and this case is known to have 
a Neel ground state which breaks spin rotation symmetry. This state has a checkerboard 
polarization of the spins, just as found in the classical ground state, and as illustrated on the 
left side of Fig. [U It can be characterized by a vector order parameter (p a which measures 
the staggered spin polarization 

V a = VtSf (1.2) 

where rji = ±1 on the two sublattices of the square lattice. In the Neel state we have 
(ip a ) 7^ 0, and we expect that the low energy excitations can be described by long wavelength 
fluctuations of a field <p a (x, r) over space, x, and imaginary time r. 

On the other hand, for A ^ 1 it is evident from Fig. [1] that the ground state preserves 
all symmetries of the Hamiltonian: it has total spin S — and can be considered to be a 
product of nearest neighbor singlet valence bonds on the J links. It is clear that this state 
cannot be smoothly connected to the Neel state, and so there must at least one quantum 
phase transition as a function A. 

Extensive quantum Monte Carlo simulations^- on this model have shown there is a direct 
phase transition between these states at a critical A c , as in Fig. [TJ The value of A c is known 
accurately, as are the critical exponents characterizing a second-order quantum phase tran- 
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sition. These critical exponents are in excellent agreement with the simplest proposal for the 
critical field theory,- which can be obtained via conventional Landau- Ginzburg arguments. 
Given the vector order parameter (p a , we write down the action in d spatial and one time 
dimension, 



as the simplest action expanded in gradients and powers of (p a which is consistent will all 
the symmetries of the lattice antiferromagnet. The transition is now tuned by varying s ~ 
(A — A c ). Notice that this model is identical to the Landau- Ginzburg theory for the thermal 
phase transition in a d + 1 dimensional ferromagnet, because time appears as just another 
dimension. As an example of the agreement: the critical exponent of the correlation length, 
u, has the same value, v = 0.711 .. ., to three significant digits in a quantum Monte Carlo 
study of the coupled dimer antiferromagnet,- and in a 5- loop analysis - of the renormalization 
group fixed point of Slg va. d — 2. Similar excellent agreement is obtained for the double- 
layer antiferromagnet^ and the coupled-plaquette antiferromagnet. 9 

In experiments, the best studied realization of the coupled-dimer antiferromagnet is 
TICUCI3. In this crystal, the dimers are coupled in all three spatial dimensions, and the 
transition from the dimerized state to the Neel state can be induced by application of pres- 
sure. Neutron scattering experiments by Ruegg and collaborators^ have clearly observed 
the transformation in the excitation spectrum across the transition, as is described by a 
simple fluctuations analysis about the mean field saddle point of Slg- in the dimerized 
phase (s > 0), a triplet of gapped excitations is observed, corresponding to the three normal 
modes of ip a oscillating about ip a = 0; as expected, this triplet gap vanishes upon approach- 
ing the quantum critical point. In a mean field analysis, valid for d > 3, the field theory 
in Eq. (jl.3p has a triplet gap of y/s. In the Neel phase, the neutron scattering detects 2 
gapless spin waves, and one gapped longitudinal mode^-. This is describe by Slg f° r s < 0, 
where (p a experiences an inverted 'Mexican hat' potential with a minumum at \(p a \ = y/\s\Jv. 
Expanding about this minimum we find that in addition to the gapless spin waves, there is 
a mode involving amplitude fluctuations of \(p a \ which has an energy gap of a/2|s|. These 
mean field predictions for the energy of the gapped modes on the two sides of the transition 
are tested in Fig. |2j the observations are in good agreement with the 1/2 exponent and the 
predicted^ \/2 ratio, providing a non-trival experimental test of the Slg field theory. 

B. Deconfined criticality 

We now consider an analog of transition discussed in Section II A\ but for a Hamiltonian 
H = H + XHi which has full square lattice symmetry at all A. For H , we choose a form 
of Hj, with Jij = J for all nearest neighbor links. Thus at A = the ground state has Neel 
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FIG. 2: Energies of the gapped collective modes across the pressure (p) tuned quantum phase 
transition in TICUCI3 observed by Ruegg et aL— . We test the description by the action Slg hi 
Eq. (|1.3p with s oc (p c — p) by comparing v2 times the energy gap for p < p c with the energy of 
the longitudinal mode for p > p c . The lines are the fits to a y/\p — p c \ dependence, testing the 1/2 
exponent. 



order, as in the left panel of Fig. [TJ We now want to choose Hi so that increasing A leads 
to a spin singlet state with spin rotation symmetry restored. A large number of choices 
have been made in the literature, and the resulting ground state invariably^ has valence 
bond solid (VBS) order; a VBS state has been observed in the organic antiferromagnet 
EtMe3P[Pd(dmit)2]2^^. The VBS state is superficially similar to the dimer singlet state 
in the right panel of Fig. [TJ the spins primarily form valence bonds with near-neighbor 
sites. However, because of the square lattice symmetry of the Hamiltonian, a columnar 
arrangement of the valence bonds as in Fig. [TJ, breaks the square lattice rotation symmetry; 
there are 4 equivalent columnar states, with the valence bond columns running along different 
directions. More generally, a VBS state is a spin singlet state, with a non-zero degeneracy 
due to a spontaneously broken lattice symmetry. Thus a direct transition between the Neel 
and VBS states involves two distinct broken symmetries: spin rotation symmetry, which is 
broken only in the Neel state, and a lattice rotation symmetry, which is broken only in the 
VBS state. The rules of Landau-Ginzburg- Wilson theory imply that there can be no generic 
second-order transition between such states. 

It has been argued that a second-order Neel- VBS transition can indeed occur—, but the 
critical theory is not expressed directly in terms of either order parameter. It involves a 
fractionalized bosonic spinor z a (a =t>4-)> an d an emergent gauge field A^. The key step is 
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to express the vector field (p a in terms of z a by 

<P a = «bZ? (1.4) 



where a a are the 2x2 Pauli matrices. Note that this mapping from (p a to z a is redundant. 
We can make a spacetime-dependent change in the phase of the z a by the field 9(x,r) 

z a -»■ e l6 z a (1.5) 

and leave ip a unchanged. All physical properties must therefore also be invariant under 
Eq. (II. 5p . and so the quantum field theory for z a has a U(l) gauge invariance, much like 
that found in quantum electrodynamics. The effective action for the z a therefore requires 
introduction of an 'emergent' U(l) gauge field (where /i = x, r is a three-component 
spacetime index). The field is unrelated the electromagnetic field, but is an internal field 
which conveniently describes the couplings between the spin excitations of the antiferromag- 
net. As we did for Slc we can write down the quantum field theory for z a and A^ by the 
constraints of symmetry and gauge invariance, which now yields 



S 7 = I d rdr 



{{dp - iA,,)z a \ +s\z a \ +u(\z a \) +—(e^xd v A 

^9 



;i.6) 



For brevity, we have now used a "relativistically" invariant notation, and scaled away the 
spin-wave velocity t>; the values of the couplings s, u are different from, but related to, those 
in Slg- The Maxwell action for A^ is generated from short distance z a fluctuations, and it 
makes A^ a dynamical field; its coupling g is unrelated to the electron charge. The action 
S z is a valid description of the Neel state for s < (the critical upper value of s will have 
fluctuation corrections away from 0), where the gauge theory enters a Higgs phase with 
(z a ) 7^ 0. This description of the Neel state as a Higgs phase has an analogy with the 
Weinberg-Salam theory of the weak interactions — in the latter case it is hypothesized that 
the condensation of a Higgs boson gives a mass to the W and Z gauge bosons, whereas here 
the condensation of z a quenches the A^ gauge boson. As written, the s > phase of S z is a 
'spin liquid' state with a S = collective gapless excitation associated with the A^ photon. 
Non-perturbative effects^ associated with the monopoles in A^ (not discussed here), show 
that this spin liquid is ultimately unstable to the appearance of VBS order. 

Numerical studies of the Neel- VBS transition have focussed on a specific lattice antifer- 
romagnet proposed by Sandvik— ~— . There is strong evidence for VBS order proximate to 
the Neel state, along with persuasive evidence of a second-order transition. However, some 
studies^ 1 ^ support a very weak first order transition. 



6 




FIG. 3: Dirac dispersion spectrum for graphene showing a 'topological' quantum phase transition 
from a hole Fermi surface for \x < to a electron Fermi surface for /i > 0. 

C. Graphene 

The last few years have seen an explosion in experimental and theoretical studies^ of 
graphene: a single hexagonal layer of carbon atoms. At the currently observed temperatures, 
there is no evident broken symmetry in the electronic excitations, and so it is not conven- 
tional to think of graphene as being in the vicinity of a quantum critical point. However, 
graphene does indeed undergo a bona fide quantum phase transition, but one without any 
order parameters or broken symmetry. This transition may be viewed as being 'topological' 
in character, and is associated with a change in nature of the Fermi surface as a function of 
carrier density. 

Pure, undoped graphene has a conical electronic dispersion spectrum at two points in the 
Brillouin zone, with the Fermi energy at the particle-hole symmetric point at the apex of the 
cone. So there is no Fermi surface, just a Fermi point, where the electronic energy vanishes, 
and pure graphene is a 'semi-metal'. By applying a bias voltage, the Fermi energy can move 
away from this symmetric point, and a circular Fermi surface develops, as illustrated in 
Fig. [3j The Fermi surface is electron-like for one sign of the bias, and hole-like for the other 
sign. This change from electron to hole character as a function of bias voltage constitutes 
the quantum phase transition in graphene. As we will see below, with regard to its dynamic 
properties near zero bias, graphene behaves in almost all respects like a canonical quantum 
critical system. 

The field theory for graphene involves fermionic degrees of freedom. Representing the 
electronic orbitals near one of the Dirac points by the two-component fermionic spinor \l/ a , 
where a is a sublattice index (we suppress spin and 'valley' indices), we have the effective 
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electronic action 

= J d 2 r J dr tft \(d T + %A T - n)5 ab + w F r x ah d x + w F r y ah d y \ ^ b 

where r* b are Pauli matrices in the sublattice space, \i is the chemical potential, v F is 
the Fermi velocity, and A T is the scalar potential mediating the Coulomb interaction with 
coupling g 2 = e 2 /e (e is a dielectric constant). This theory undergoes a quantum phase 
transition as a function of fi, at fi = 0, similar in many ways to that of Slg as a function 
of s. The interaction between the fermionic excitations here has coupling g 2 , which is the 
analog of the non-linearity u in Slg- However, while u scaled to a fixed point non-zero fixed 
point value under the renormalization group flow, g flows logarithmically slowly to zero. For 
many purposes, it is safe to ignore this flow, and to set g equal to fixed value; this value is 
characterized by the dimensionless 'fine structure constant' a = g 2 /(hvp) which is of order 
unity in graphene. 



D. Spin density waves 

Finally, we consider the onset of Neel order, as in Section II A\ but in a metal rather 
than an insulator. It is conventional to refer to such metallic Neel states as have spin 
density wave (SDW) order. Our discussion here is motivated by application to the cuprate 
superconductors: there is good evidence^ 1 ^ that the transition we describe below is present 
in the electron-doped cuprates, and proposals of its application to the hole-doped cuprates 
will be discussed in Section ED 

We begin with the band structure describing the cuprates in the over-doped region, well 
away from the Mott insulator. Here the electrons q q are described by the Hamiltonian 

H c — ~ ^ j tijC ia Ci a = ^ ] £k c ka C ka 
i<j k 

with tij non-zero for first, second and third neighbors on the square lattice. This leads to 
Fermi surface shown in the right-most panel of Fig. HI As in Section II CI we will consider 
topological changes to this Fermi surface, but is induced here by a conventional SDW (Neel) 
order parameter ip a . From the structure of the ordering in Fig. [1] we see that (p a carries 
momentum K = (jr, ir), and so will transfer momentum K while scattering the c a fermions; 
this leads to the coupling 

#sdw = V a C k,a CT a/3 C k+K,/3- 
k,a,/3 

We can now follow the evolution of the Fermi surface under the onset of SDW order by 
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^ Increasing SDW order 

^ 

FIG. 4: Evolution of the Fermi surface of the hole doped cuprates in a conventional SDW theory^ 
as a function of the magnitude of the SDW order | ip a \ . The right panel is the large Fermi surface 
state with no SDW order, with states contiguous to k = occupied by electrons. The onset of 
SDW order induces the formation of electron (red) and hole (blue) pockets. With further increase 
of \(p a \, the electron pockets disappear and only hole pockets remain (the converse happens in the 
last step for the electron-doped cuprates. 

diagonalizing H c + H s< i w for constant (p a : the results of this^ are shown in Fig. HI The 
second panel from the right shows the Fermi surface obtained by translating the original 
Fermi surface by K, and the remaining panels show the consequences of mixing between the 
states at momentum k and k + K. 

We are now interested in the nature of the quantum critical point at the onset of SDW 
order, indicated by the vertical line in Fig. HJ Note that this transition combines the features 
of Sections II Al and ll CI it has an order parameter, as in the coupled dimer antiferromagnet, 
and it has a topological change in the Fermi surface, as in graphene. Thus a complete theory 
is just the combination of these two features: Slg + H c + H s <i w . A common assert ion2&>2£ 
is that we can pay less attention to the Fermi surface change by simply integrating out the 
c a fermions and working with the resulting modified action for (f a . Right at the quantum 
critical point, the SDW fluctuations, ip a connect points on the large Fermi surface, and so 
can decay into a large density of states of particle-hole excitations. The damping induced 
by this particle-hole continuum modifies the effective action for (p a from Eq. (jl.3p by adding 
a relevant dissipative term:— >2I 

Sn = S LG + J^f ^M|^(*,a>)l 2 (1.8) 

More recent analyses 38 have indicated that this procedure is likely not correct in two spatial 
dimensions, and that the Fermi surface change has to be treated in a more fundamental 
manner. 
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FIG. 5: Finite temperature crossovers of the coupled dimer antiferromagnet in Fig. [TJ 



II. FINITE TEMPERATURE CROSSOVERS 



The previous section has described four model systems at T = 0: we examined the change 
in the nature of the ground state as a function of some tuning parameter, and motivated 
a quantum field theory which describes the low energy excitations on both sides of the 
quantum critical point. 

We now turn to the important question of the physics at non-zero temperatures. All of 
the models share some common features, which we will first explore for the coupled dimer 
antiferromagnet. For A > A c (or s > in Slg), the excitations consist of a triplet of S — 1 
particles (the 'triplons'), which can be understood perturbatively in the large A expansion 
as an excited S = 1 state on a dimer, hopping between dimers (see Fig. [5]). The mean field 
theory tells us that the excitation energy of this dimer vanishes as \fs upon approaching 
the quantum critical point. Fluctuations beyond mean field, described by Slg, show that 
the exponent is modified to s zu , where z — 1 is the dynamic critical exponent, and v is 
the correlation length exponent. Now imagine turning on a non-zero temperature. As long 
as T is smaller than the triplon gap, i.e. T < s zu , we expect a description in terms of a 
dilute gas of thermally excited triplon particles. This leads to the behavior shown on the 
right-hand-side of Fig. [51 delimited by the crossover indicted by the dashed line. Note that 
the crossover line approaches T = only at the quantum critical point. 

Now let us look a the complementary behavior at T > on the Neel-ordered side of the 
transition, with s < 0. In two spatial dimensions, thermal fluctuations prohibit the breaking 
of a non-Abelian symmetry at all T > 0, and so spin rotation symmetry is immediately 
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restored. Nevertheless, there is an exponentially large spin correlation length, £, and at 
distances shorter than £ we can use the ordered ground state to understand the nature of 
the excitations. Along with the spin-waves, we also found the longitudinal 'Higgs' mode 
with energy v '—2s in mean field theory. Thus, just as was this case for s > 0, we expect this 
spin-wave+Higgs picture to apply at all temperatures lower than the natural energy scale; 
i.e. for T < {—s) zu . This leads to the crossover boundary shown on the left-hand-side of 
Fig. EJ 

Having delineated the physics on the two sides of the transition, we are left with the 
crucial quantum critical region in the center of Fig. [51 This is present for T > \s\ zu , i.e. at 
higher temperatures in the vicinity of the quantum critical point. To the left of the quantum 
critical region, we have a description of the dynamics and transport in terms of an effectively 
classical model of spin waves: this is the 'renormalized classical' regime of Ref.— . To the 
right of the quantum critical region, we again have a regime of classical dynamics, but now in 
terms of a Boltzmann equation for the triplon particles. A key property of quantum critical 
region is that there is no description in terms of either classical particles or classical waves at 
the times of order the typical relaxation time, r r , of thermal excitations. Instead, quantum 
and thermal effects are equally important, and involve the non-trivial dynamics of the fixed- 
point theory describing the quantum critical point. Note that while the fixed-point theory 
applies only at a single point (A = A c ) at T = 0, its influence broadens into the quantum 
critical region at T > 0. Because there is no characteristing energy scale associated with 
the fixed-point theory, k%T is the only energy scale available to determine r r at non-zero 
temperatures. Thus, in the quantum critical region 

k B T 

where C is a universal constant dependent only upon the universality class of the fixed point 
theory i.e. it is universal number just like the critical exponents. This value of r r determines 
the 'friction coefficients' associated with the dissipative relaxation of spin fluctuations in the 
quantum critical region. It is also important for the transport co-efficients associated with 
conserved quantities, and this will be discussed in Section IHIl 

Let us now consider the similar T > crossovers for the other models of Section [H 
The Neel-VBS transition of Section II Bl has crossovers very similar to those in Fig. 
with one important difference. The VBS state breaks a discrete lattice symmetry, and 
this symmetry remains broken for a finite range of non-zero temperatures. Thus, within 
the right-hand 'triplon gas' regime of Fig. [5j there is a phase transition line at a critical 
temperature Tvbs- The value of Tvbs vanishes very rapidly as s \ 0, and is controlled by 
the non-perturbative monopole effects which were briefly noted in Section II Bl 

For graphene, the discussion above applied to Fig. [3] leads to the crossover diagram shown 
in Fig. [61 as noted by Sheehy and Schmalian^i. We have the Fermi liquid regimes of the 
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FIG. 6: Finite temperature crossovers of graphene as a function of electron density n (which is 



tuned by \x in Eq. (|1.7p ) and temperature, T. Adapted from Ref. 131 



electron- and hole- like Fermi surfaces on either side of the critical point, along with an 
intermediate quantum critical Dirac liquid. A new feature here is related to the logarithmic 
flow of the dimensionless 'fine structure constant' a controlling the Coulomb interactions, 
which was noted in Section II CI In the quantum critical region, this constant takes the 
typical value a ~ l/ln(l/T). Consequently for the relaxation time in Eq. (12. ip we have 
C ~ ln 2 (l/T). This time determines both the width of the electron spectral functions, and 
also the transport co-efficients, as we will see in Section IIHI 

Finally, we turn to the spin density wave transition of Section II Dl From the evolution 
in Fig. H] and the discussion above, we have the crossover phase diagram™ in Fig. [7J As in 
Fig. El there is no transition at non-zero temperatures, and now the crossover is between the 
topologically distinct Fermi surface configurations of Fig. |U The Hertz action for the SDW 
fluctuations in Eq. fll.Sjl predicts^ 1 ^ logarithmic corrections to the leading scaling behavior, 
similar to those for graphene. The interplay of such SDW fluctuations with the topological 
change in the Fermi surface configuration is not fully understood^, and is labeled as the 
'strange metal' regime in Fig. [7J 



III. QUANTUM CRITICAL TRANSPORT 

We now turn to the 'transport' properties in the quantum critical region: we consider the 
response functions associated with any globally conserved quantity. For the antiferromag- 
netic systems in Sections II Al and II B\ this requires consideration of the transport of total 
spin, and the associated spin conductivities and diffusivities. For graphene, we can consider 
charge and momentum transport. Our discussion below will also apply to the superfluid- 
insulator transition: for bosons in a periodic potential, this transition is described^ by a field 
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FIG. 7: Finite temperature crossovers of near the SDW ordering transition of Fig. UJ Here, looking 
to application to the cuprates to be discussed in Section IVIl we have assumed that the carrier 
density, x, tunes the system across the SDW transition. 



theory closely related to that in Eq. (11. 3p . However, we will primarily use a language appro- 
priate to charge transport in graphene below. We will describe the properties of a generic 
strongly-coupled quantum critical point and mention, where appropriate, the changes due 
to the logarithmic flow of the coupling in graphene. 

In traditional condensed matter physics, transport is described by identifying the low- 
lying excitations of the quantum ground state, and writing down 'transport equations' for 
the conserved charges carried by them. Often, these excitations have a particle-like nature, 
such as the 'triplon' particles of Fig. [5] or the electron or hole quasiparticles of the Fermi 
liquids in Fig. [6j In other cases, the low-lying excitations are waves, such as the spin- 
waves in Fig. [5j and their transport is described by a non-linear wave equation (such as the 
Gross-Pitaevski equation). However, as we have discussed in Section [III neither description 
is possible in the quantum critical region, because the excitations do not have a particle-like 
or wave-like character. 

Despite the absence of an intuitive description of the quantum critical dynamics, we can 
expect that the transport properties should have a universal character determined by the 
quantum field theory of the quantum critical point. In addition to describing single excita- 
tions, this field theory also determines the S"-matrix of these excitations by the renormal- 
ization group fixed-point value of the couplings, and these should be sufficient to determine 
transport properties 35 . The transport co-efficients, and the relaxation time to local equilib- 
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rium, are not proportional to a mean free scattering time between the excitations, as is the 
case in the Boltzmann theory of quasiparticles. Such a time would typically depend upon 
the interaction strength between the particles. Rather, the system behaves like a "perfect 
fluid" in which the relaxation time is as short as possible, and is determined universally 
by the absolute temperature, as indicated in Eq. (12. ip . Indeed, it was conjectured in Ref.- 
that the relaxation time in Eq. (12.1 p is a generic lower bound for interacting quantum sys- 
tems. Thus the non-quantum-critical regimes of all the phase diagrams in Section [TT] have 
relaxation times which are all longer than Eq. (12. ip . 

The transport co-efficients of this quantum-critical perfect fluid also do not depend upon 
the interaction strength, and can be connected to the fundamental constants of nature. In 
particular, the electrical conductivity, a, is given by (in two spatial dimensions)^ 5 - 

<t q = (3.1) 

where $ CT is a universal dimensionless constant of order unity, and we have added the sub- 
script Q to emphasize that this is the conductivity for the case of graphene with the Fermi 
level at the Dirac point (for the superfluid-insulator transition, this would correspond to 
bosons at integer filling) with no impurity scattering, and at zero magnetic field. Here e* 
is the charge of the carriers: for a superfluid-insulator transition of Cooper pairs, we have 
e* = 2e, while for graphene we have e* = e. The renormalization group flow of the 'fine struc- 
ture constant' a of graphene to zero at asymptotically low T, allows an exact computation in 
this case^: $ CT ~ 0.05 ln 2 (l/T). For the superfluid-insulator transition, <$> a is T-independent 
(this is the generic situation with non-zero fixed point values of the interactional) but it has 
only been computed 1,35 to leading order in expansions in 1/N (where N is the number of 
order parameter components) and in 3 — d (where d is the spatial dimensionality). However, 
both expansions are neither straightforward nor rigorous, and require a physically motivated 
resummation of the bare perturbative expansion to all orders. It would therefore be valuable 
to have exact solutions of quantum critical transport where the above results can be tested, 
and we turn to such solutions in the next section. 

In addition to charge transport, we can also consider momentum transport. This was 
considered in the context of applications to the quark-gluon plasma 3 ^; application of the 
analysis of Ref.— shows that the viscosity, r], is given by 

V ~ = (3-2) 

s k B 

where s is the entropy density, and again § v is a universal constant of order unity. The value 
of ^ has recently been computed^- for graphene, and again has a logarithmic T dependence 
because of the marginally irrelevant interaction: <& v w 0.008 ln 2 (l/T). 

We conclude this section by discussing some subtle aspects of the physics behind the 
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seemingly simple result quantum-critical in Eq. (13. ip . For simplicity, we will consider the 
case of a "relativistically" invariant quantum critical point in 2+1 dimensions (such as the 
field theories of Section II Al and II B[ but marginally violated by graphene, a subtlety we 
ignore below). Consider the retarded correlation function of the charge density, x(k,ou), 
where k = |k| is the wavevector, and u is frequency; the dynamic conductivity, cr(u), is 
related to x by the Kubo formula, 

a(u) = lim -i-x(*> w). (3.3) 

It was argued in Ref.— that despite the absence of particle-like excitations of the critical 
ground state, the central characteristic of the transport is a crossover from collisionless to 
collision-dominated transport. At high frequencies or low temperatures, the limiting form 
for x reduces to that at T = 0, which is completely determined by relativistic and scale 
invariance and current conversion upto an overall constant 

*2 / 2 *2 

X {k : u) = e —K : , a{u) = e —K ; fa » k B T, (3.4) 

11 y V k l — [cj + ir\Y "> 

where K is a universal number—. However, phase-randomizing collisions are intrinsically 
present in any strongly interacting critical point (above one spatial dimension) and these 
lead to relaxation of perturbations to local equilibrium and the consequent emergence of 
hydrodynamic behavior. So at low frequencies, we have instead an Einstein relation which 
determines the conductivity with 

Dk 2 e* 2 

x (k,u) = e* 2 Xc T r r2 , a(uj) = e* 2 Xc D = —Q^i ; fa<^k B T, (3.5) 

Uk z — too n 

where Xc is the compressibility and D is the charge diffusion constant. Quantum critical 
scaling arguments show that the latter quantities obey 

h^v z k B l 

where 0^2 are universal numbers. A large number of papers in the literature, particularly 
those on critical points in quantum Hall systems, have used the collisionless method of 
Eq. (13 .4p to compute the conductivity. However, the correct d.c. limit is given by Eq. (13. 5p . 
and the universal constant in Eq. (13. ip is given by $ CT = ©102- Given the distinct physical 
interpretation of the collisionless and collision-dominated regimes, we expect that K ^ @i@2- 
This has been shown in a resummed perturbation expansion for a number of quantum critical 
points 1 . 
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IV. EXACT RESULTS FOR QUANTUM CRITICAL TRANSPORT 



The results of Section UTT1 were obtained by using physical arguments to motivate resum- 
mations of perturbative expansions. Here we shall support the ad hoc assumptions behind 
these results by examining an exactly solvable model of quantum critical transport. 

The solvable model may be viewed as a generalization of the gauge theory in Eq. (11. 6p 
to the maximal possible supersymmetry. In 2+1 dimensions, this is known as M = 8 
supersymmetry. Such a theory with the U(l) gauge group is free, and so we consider the 
non-Abelian Yang-Millis theory with a SU(iV) gauge group. The resulting supersymmetric 
Yang-Mills (SYM) theory has only one coupling constant, which is the analog of the electric 
charge g in Eq. (jl.fip . The matter content is naturally more complicated than the complex 
scalar z a in Eq. (II. 6p . and also involves relativistic Dirac fermions as in Eq. ( II. 7p . However 
all the terms in the action for the matter fields are also uniquely fixed by the single coupling 
constant g. Under the renormalization group, it is believed that g flows to an attractive 
fixed point at a non-zero coupling g = g*; the fixed point then defines a supersymmetric 
conformal field theory in 2+1 dimensions (a SCFT3), and we are interested here in the 
transport properties of this SCFT3. 

A remarkable recent advance has been the exact solution of this SCFT3 in the N — > oo 
limit using the AdS/CFT correspondence^. The solution proceeds by a dual formulation 
as a four- dimensional supergravity theory on a spacetime with uniform negative curvature: 
anti-de Sitter space, or AdS 4 . Remarkably, the solution is also easily extended to non-zero 
temperatures, and allows direct computation of the correlators of conserved charges in real 
time. At T > a black hole appears in the gravity, resulting in an AdS-Schwarzschild 
spacetime, and T is also the Hawking temperature of the black hole; the real time solutions 
also extend to T > 0. 

The results of a full computation^ of the density correlation function, x(k, are shown 
in Fig. [S] and O The most important feature of these results is that the expected limiting 
forms in the collisionless (Eq. (13. 4p ) and collision-dominated (Eq. (13. 5p are obeyed. Thus the 
results do display the collisionless to collision-dominated crossover at a frequency of order 
ksT/h, as was postulated in Section HTT1 

An additional important feature of the solution is apparent upon describing the full 
structure of both the density and current correlations. Using spacetime indices (/i, v — t,x, y) 
we can represent these as the tensor XiwO*-, oS), where the previously considered x = Xtt- At 
T > 0, we do not expect x^u to be relativistically covariant, and so can only constrain it by 
spatial isotropy and density conservation. Introducing a spacetime momentum p^ = (uj, k), 
and setting the velocity v = 1, these two constraints lead to the most general form 

X>(k,w) = e —^(p T ^K T (k } u) + P^K L (k^)) (4.1) 
where p 2 = rf v p^p v with T]^ = diag(— 1, 1,1), and P^ u and Ph, are orthogonal projectors 
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FIG. 8: Spectral weight of the density correlation function of the SCFT3 with N = 8 supersym- 
metry in the collisionless regime. 
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FIG. 9: As in Fig. [8l but for the collision-dominated regime. 



defined by 

pT _ pT _ pT _ q p T _ x. _ Mi P L - (n - 2n^L\ _ / 4 9) 

MM — MM iO ' My ~~ °y p ' A 1 ^ ~~ ylV ^2 J ysw> M"M 

with the indices i,j running over the 2 spatial components. The two functions K T,L (k,u) 
define all the correlators of the density and the current, and the results in Eqs. (13.51) and 
(13. 4p are obtained by taking suitable limits of these functions. We will also need below the 
general identity 

K T (0,oo) = K L (0,oo), (4.3) 
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which follows from the analyticity of the T > current correlations at k = 0. 

The relations of the previous paragraph are completely general and apply to any theory. 
Specializing to the AdS-Schwarzschild solution of SYM3, the results were found to obey a 
simple and remarkable identity 43 : 

K L (k,u)K T (k,u) = K? (4.4) 

where /C is a known pure number, independent of u and k. It was also shown that such 
a relation applies to any theory which is equated to classical gravity on AdS4, and is a 
consequence of the electromagnetic self-duality of its four-dimensional Maxwell sector. The 
combination of Eqs. (14.31) and (14.41) fully determines the Xia> correlators at k = 0: we find 
K L (0,u) = K t (0,uj) = /C, from which it follows that the k = conductivity is frequency 
independent and that $ CT = O1O2 — K — JC. These last features are believed to be special 
to theories which are equivalent to classical gravity, and not hold more generally. 

We can obtain further insight into the interpretation of Eq. (14. 4p by considering the 
field theory of the superfluid-insulator transition of lattice bosons at integer filling. As 
we noted earliear, this is given by the field theory in Eq. (II. 3p with the field (p a having 2 
components. It is known that this 2-component theory of relativistic bosons is equivalent 
to a dual relativistic theory, S of vortices, under the well-known 'particle-vortex' duality^ 4 -. 
Ref.— considered the action of this particle-vortex duality on the correlation functions in 
Eq. (14.11) . and found the following interesting relations: 

K L (k,ou)K T (k,Lo) = 1 , K T (k,ou)K L (k,Lo) = 1 (4.5) 

where K L ' T determine the vortex current correlations in S as in Eq. (14.11) . Unlike Eq. (14.41) . 
Eq. ( 14.51) does noHully determine the correlation functions at k = 0: it only serves to reduce 
the 4 unknown functions K L,T , K L,T to 2 unknown functions. The key property here is that 
while the theories Slg an d S are dual to each other, they are not equivalent, and the theory 
Slg is n °t self-dual. 

We now see that Eq. (14. 4p implies that the classical gravity theory of SYM3 is self-dual 
under an analog of particle- vortex duality 43 . It is not expected that this self-duality will 
hold when quantum gravity corrections are included; equivalently, the SYM3 at finite N 
is expected to have a frequency dependence in its conductivity at k = 0. If we apply the 
AdS/CFT correspondence to the superfluid-insulator transition, and approximate the latter 
theory by classical gravity on AdS4, we immediately obtain the self-dual prediction for the 
conductivity, $ CT = 1. This value is not far from that observed in numerous experiments, 
and we propose here that the AdS / CFT correspondence offers a rationale for understanding 
such observations. 
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V. HYDRODYNAMIC THEORY 



The successful comparison between the general considerations of Section IHIt and the 
exact solution using the AdS/CFT correspondence in Section HVl emboldens us to seek a 
more general theory of low frequency (fiw <C ksT) transport in the quantum critical regime. 
We will again present our results for the special case of a relativistic quantum critical point in 
2+1 dimensions (a CFT3), but it is clear that similar considerations apply to a wider class of 
systems. Thus we can envisage applications to the superfluid-insulator transition, and have 
presented scenarios under which such a framework can be used to interpret measurements of 
the Nernst effect in the cuprates^ -. We have also described a separate set of applications to 
graphene 3 ^: while graphene is strictly not a CFT3, the Dirac spectrum of electrons leads to 
many similar results, especially in the inelastic collision-dominated regime associated with 
the quantum critical region. These results on graphene are reviewed in a separate paper—, 
where explicit microscopic computations are also discussed. 

Our idea is to relax the restricted set of conditions under which the results of Section II III 
were obtained. We will work within the quantum critical regimes of the phase diagrams of 
Section [Til but now allow a variety of additional perturbations. First, we will move away 
from the particle-hole symmetric case, allow a finite density of carriers. For graphene, 
this means that \x is no longer pinned at zero; for the antiferromagnets, we can apply an 
external magnetic field; for the superfluid-insulator transition, the number density need not 
be commensurate with the underlying lattice. For charged systems, such as the superfluid- 
insulator transition or graphene, we allow application of an external magnetic field. Finally, 
we also allow a small density of impurities which can act as a sink of the conserved total 
momentum of the CFT3. In all cases, the energy scale associated with these perturbations 
is assumed to be smaller than the dominant energy scale of the quantum critical region, 
which is ksT. The results presented below were obtained in two separate computations, 
associated with the methods described in Sections II I II and IIV} and are described in the two 
subsections below. 

A. Relativistic magnetohydrodynamics 

With the picture of relaxation to local equilibrium at frequencies ftw <C k B T developed 
in Ref. 35 , we postulate that the equations of relativistic magnetohydrodynamics should de- 
scribe the low frequency transport. The basic principles involved in such a hydrodynamic 
computation go back to the nineteenth century: conservation of energy, momentum, and 
charge, and the constraint of the positivity of entropy production. Nevertheless, the re- 
quired results were not obtained until our recent work 40 : the general case of a CFT3 in 
the presence of a chemical potential, magnetic field, and small density of impurities is very 
intricate, and the guidance provided by the dual gravity formulation was very helpful to us. 
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In this approach, we do not have quantitative knowledge of a few transport co-efficients, and 
this is complementary to our ignorance of the effective couplings in the dual gravity theory 
to be discussed in Section IV Bl 

The complete hydrodynamic analysis can be found in Ref.— . The analysis is intricate, 
but is mainly a straightforward adaption of the classic procedure outlined by Kadanoff and 
Martin^ to the relativistic field theories which describe quantum critical points. We list the 
steps: 

1. Identify the conserved quantities, which are the energy-momentum tensor, T^, and 
the particle number current, J M . 

2. Obtain the real time equations of motion, which express the conservation laws: 

d u T^ = F^J u , 9^ = 0; (5.1) 

here F^ u represents the externally applied electric and magnetic fields which can 
change the net momentum or energy of the system, and we have not written a term 
describing momentum relaxation by impurities. 

3. Identify the state variables which characterize the local thermodynamic state — we 
choose these to be the density, p, the temperature T, and an average velocity u^. 

4. Express T^ u and J 71 in terms of the state variables and their spatial and temporal 
gradients; here we use the properties of the observables under a boost by the velocity 
u M , and thermodynamic quantities like the energy density, e, and the pressure, P, 
which are determined from T and p by the equation of state of the CFT. We also 
introduce transport co-efficients associated with the gradient terms. 

5. Express the equations of motion in terms of the state variables, and ensure that the 
entropy production rate is positive^. This is a key step which ensures relaxation to 
local equilibrium, and leads to important constraints on the transport co-efficients. In 
d = 2, it was found that situations with the velocity spacetime independent are 
characterized by only a single independent transport co-efficient^. This we choose to 
be the longitudinal conductivity at B = 0. 

6. Solve the initial value problem for the state variables using the linearized equations of 
motion. 

7. Finally, translate this solution to the linear response functions, as described in Ref.— . 
B. Dyonic black hole 

Given the success of the AdS / CFT correspondence for the specific supersymmetric model 
in Section HV] we boldly assume a similar correspondence for a generic CFT3. We assume 
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that each CFT3 is dual to a strongly-coupled theory of gravity on AdS4. Furthermore, 
given the operators associated with the perturbations away from the pure CFT3 we want 
to study, we can also deduce the corresponding perturbations away from the dual gravity 
theory. So far, this correspondence is purely formal and not of much practical use to us. 
However, we now restrict our attention to the hydrodynamic, collision dominated regime, 
hu <C k B T of the CFT3. We would like to know the corresponding low energy effective 
theory describing the quantum gravity theory on AdS4. Here, we make the simplest possible 
assumption: the effective theory is just the Einstein-Maxwell theory of general relativity and 
electromagnetism on AdS 4 . As in Section HV] the temperature T of CFT3 corresponds to 
introducing a black hole on AdS 4 whose Hawking temperature is T. The chemical potential, 
fi, of the CFT3 corresponds to an electric charge on the black hole, and the applied magnetic 
field maps to a magnetic charge on the black hole. Such a dynoic black hole solution of the 
Einstein-Maxwell equations is, in fact, known: it is the Reissner-Nordstrom black hole. 

We solved the classical Einstein-Maxwell equations for linearized fluctuations about the 
metric of a dyonic black hole in a space which is asymptotically AdS4. The results were used 
to obtain correlators of a CFT3 using the prescriptions of the AdS/CFT mapping. As we 
have noted, we have no detailed knowledge of the strongly-coupled quantum gravity theory 
which is dual to the CFT3 describing the superfluid-insulator transition in condensed matter 
systems, or of graphene. Nevertheless, given our postulate that its low energy effective field 
theory essentially captured by the Einstein-Maxwell theory, we can then obtain a powerful 
set of results for CFT3s. 



C. Results 

In the end, we obtained complete agreement between the two independent computations 
in Sections IV Al and IV B\ after allowing for their distinct equations of state. This agree- 
ment demonstrates that the assumption of a low energy Einstein-Maxwell effective field 
theory for a strongly coupled theory of quantum gravity is equivalent to the assumption of 
hydrodynamic transport for hu <C k^T in a strongly coupled CFT3. 

Finally, we turn to our explicit results for quantum critical transport with hu <C k%T ' . 

First, consider adding a chemical potential, /i, to the CFT3. This will induce a non-zero 
number density of carriers p. The value of p is defined so that the total charge density 
associated with p is e*p. Then the electrical conductivity at a frequency u is 

e* 2 e*Vf 2 1 

a(u) = —$a + 7 n w ; r (5.2) 

1 ; h " (e + P) {-ico + l/r imp ) 1 ; 

In this section, we are again using the symbol v to denote the characteristic velocity of the 
CFT3 because we will need c for the physical velocity of light below. Here e is the energy 
density and P is the pressure of the CFT3. We have assumed a small density of impurities 
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which lead to a momentum relaxation time T ira ^r^-. In general, p, e, P, and l/r imp 
will be functions of pjk^T which cannot be computed by hydrodynamic considerations 
alone. However, apart from these quantities are usually amenable to direct perturbative 
computations in the CFT3, or by quantum Monte Carlo studies. The physical interpretation 
of Eq. (15. 2p should be evident: adding a charge density p leads to an additional Drude-like 
contribution to the conductivity. This extra current cannot be relaxed by collisions between 
the unequal density of particle and hole excitations, and so requires an impurity relaxation 
mechanism to yield a finite conductivity in the d.c. limit. 

Now consider thermal transport in a CFT3 with a non-zero p. The d.c. thermal conduc- 
tivity, k, is given by 

mwi 

in the absence of impurity scattering, l/r imp — > 0. This is a Wiedemann- Franz-like relation, 
connecting the thermal conductivity to the electrical conductivity in the p = CFT. Note 
that k diverges as p — > 0, and so the thermal conductivity of the p = CFT is infinite. 

Next, turn on a small magnetic field B; we assume that B is small enough that the spacing 
between the Landau levels is not as large as ksT. The case of large Landau level spacing is 
also experimentally important, but cannot be addressed by the present analysis. Initially, 
consider the case p = 0. In this case, the result Eq. (I5.3P for the thermal conductivity is 
replaced by 

K = *7 (~/r) {k B TB/(hc/e*)) (5 - 4) 

also in the absence of impurity scattering, l/r imp — > 0. This result relates k to the electrical 
resistance at criticality, and so can be viewed as Wiedemann-Franz-like relation for the 
vortices. A similar 1/B 2 dependence of k appeared in the Boltzmann equation analysis of 
Ref.—, but our more general analysis applies in a wider and distinct regime 36 , and relates 
the co-efficient to other observables. 

We have obtained a full set of results for the frequency- dependent thermo-electric response 
functions at non-zero B and p. The results are lengthy and we refer the reader to Ref.— for 
explicit expressions. Here we only note that the characteristic feature^^ of these results 
is a new hydrodynamic cyclotron resonance. The usual cyclotron resonance occurs at the 
classical cyclotron frequency, which is independent of the particle density and temperature; 
further, in a Galilean-invariant system this resonance is not broadened by electron-electron 
interactions alone, and requires impurities for non-zero damping. The situtation for our 
hydrodynamic resonance is very different. It occurs in a collision-dominated regime, and its 
frequency depends on the density and temperature: the explicit expression for the resonance 
frequency is 



22 



Further, the cyclotron resonance involves particle and hole excitations moving in opposite 
directions, and collisions between them can damp the resonance even in the absence of 
impurities. Our expression for this intrinsic damping frequency 



relating it to the quantum-critical conductivity as a measure of collisions between counter- 
propagating particles and holes. We refer the reader to a separate discussion 36 of the exper- 
imental conditions under which this hydrodynamic cyclotron resonance may be observed. 

VI. THE CUPRATE SUPERCONDUCTORS 

We close this article by mentioning application to the cuprate superconductors. The 
phenomenology of these materials is very involved, and they have so far resisted a global in- 
terpretation in terms of the canonical quantum-critical crossover phase diagrams discussed 
here. Surely, one of the important complexities is that they involve at least two order param- 
eters: those associated with c?-wave superconductivity and with SDWs (both commensurate 
and incommensurate with the underlying lattice). In addition there are also topological 
changes in the Fermi surface discussed in Section II Dl and Fig. HJ 

Here we mention recent ideas^~— which attempt a synthesis based upon a combination 
of the crossover phase diagrams we have discussed here. The ideas build upon the results 
of a number of recent experiments^ - — which have explored the phase diagram of both the 
hole- and electron-doped cuprates in a magnetic field applied perpendicular to the layers. 
We do not wish to enter into the intricate details of the structure of the phase diagram in 
the field-doping plane, and refer the reader to recent reviews^ 3 - 1 ^. We will limit ourselves to 
the assertion that the experiments support the proposal^ that the underlying metallic state 
of the cuprates undergoes an SDW transition as in Fig. [7J along with the topological change 
of the Fermi surface, which was illustrated in Fig. H] for the simplest case of commensurate 
ordering. Recent computations^ and observations^ of thermoelectric effects have provided 
additional support for this proposal. 

The main ingredient missing from Fig.[7Jis c?-wave superconductivity. It is known that the 
SDW fluctuations can induce <i-wave pairing for the Fermi surface configuration shown in 
Fig. HI So let us superimpose a superconducting region on Fig. [7J and this leads us to Fig. [TUl 
Arguments that this superconducting region will have the dome shape around the quantum 
critical point were made in Ref.— . Here we wish to draw attention to a crucial new feature 
shown in Fig. (TUJ The onset of superconductivity shifts the spin ordering quantum critical 
point from its apparent position in the underlying metal at x = x m , to its actual position in 
the superconductor at x = x s . This shift is indicated by the large arrow in Fig. [10J We have 
argued^ - — that this shift helps resolve a number of long-standing puzzles in understanding 
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FIG. 10: Influence of the onset of superconductivity on the crossovers near the SDW ordering 
transition in Fig. [7J 



the experimental observations in the cuprates in terms of the quantum-critical crossover 
phase diagrams. 

What is the physical origin of the shift in the position of the quantum critical point? In 
phenomenological terms, this shift can be understood as a competition between the SDW and 
superconducting order parameters^. Once there is an onset of superconductivity, it repels 
the incipient SDW ordering, and shifts it away from the superconducting region to lower 
carrier concentration. A microscopic theory has also been presented in Ref.~: it requires 
special attention to the physics of the Fermi surface, and the fact that SDW ordering and 
superconductivity are 'eating up' the same portions of the Fermi surface. Such phenomena 
are especially pronounced near the Mott insulating state at zero doping. 

Once the shift has occurred, the crossovers near the spin-ordering transition within the 
superconductor are especially simple. There are no longer any Fermi surfaces to contend 
with, and so we can largely focus on the fluctuations of the order parameter. The latter 
reduce^ 3 to those presented in Fig. [51 with critical theories like those presented in Sections II Al 
and HI 

So the bottom line is that the proposed cuprate phase diagram in Fig. [10] is a combination 
of the crossovers in Fig. and [7J along with the all-important shift in the position of the 
quantum critical point. 
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It is now interesting to examine the subtle evolution in the physics as the temperature is 
lowered for x s < x < x m . Above T c , we have a metallic state with fluctuating SDW order 
whose strength increases as T is lowered. Once we cross T c , the growth in SDW ordering is 
arrested (because of the competition with superconductivity); the system is now to the right 
of the SDW ordering quantum critical point, and starts to recover aspects of the physics of 
the large Fermi surface. However, it has been argued 63 * that forms of translational symmetry 
breaking, such as VBS or nematic ordering may appear as T — > for x s < x < x m . In this 
case, the left crossover line emerging from x m in Fig. [10] would become an actual thermal 
phase transition at which lattice symmetry is restored. 
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